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Abstract 

We found that three types of tethered surface model undergo a first-order phase 
transition between the smooth and the crumpled phase. The first and the third are 
discrete models of Helfrich, Polyakov, and Kleinert, and the second is that of Nambu 
and Goto. These are curvature models for biological membranes including artificial 
vesicles. The results obtained in this paper indicate that the first-order phase tran- 
sition is universal in the sense that the order of the transition is independent of 
discretization of the Hamiltonian for the tethered surface model. 
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1 Introduction 



A considerable number of attempts have been made at clarifying the phase 
structure of the surface model of Helfrich, Polyakov, and Kleinert [1,2,3,4,5,6,7,8,9,10,11,12,13,14] 
henceforce we will refer to this as the HPK model. Tethered surface models are 
defined on triangulated fixed connectivity surfaces representing polymerized 
biological membranes or membranes in the gel phase [7], and they are classi- 
fied into a major class of the HPK model [15,16,17,18,19,20,21,22,23,24,25,26]. 
Fluid surface models are considered a different class of the HPK model defined 
on dynamically triangulated surfaces representing these biological membranes 
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in the fluid phase, however, we will not discuss the fluid surface model in this 
paper [27,28,29,30,31,32,33,34]. 

An interesting problem that remains unanswered at present is whether the 
tethered self-intersecting (or phantom) surface model undergoes a continuous 
transition or a discontinuous one. The model is expected to undergo a phase 
transition separating the smooth phase in the limit b-^oo from the crumpled 
one in the limit 6 — > 0. The mean field analysis predicts that the model un- 
dergoes a discontinuous phase transition [12], whereas the large-D expansion 
predicts a continuous transition [11]. Biological membranes including artifi- 
cial vesicles were thought to undergo the crumpling transition. However, we 
currently have no experimental evidence for such a transition in biological 
membranes, except for an investigation on artificial vesicles [35]. 

In a recent numerical study, it was reported that the tethered model of HPK 
undergoes a first-order phase transition on spherical surfaces relatively larger 
than those used so far [26] . The bending energy and the Gaussian term used in 
[26] are those widely accepted as standard discrete Hamiltonians in numerical 
studies on the HPK model. The bending energy is of the form 1— rij • , where 
iij is the unit normal vector of the triangle i. The Gaussian term is given by 
the sum of bond length squares. 

However, it has not yet been clarified whether the order of phase transition 
depends on the discretization of the Hamiltonian. If the order of phase transi- 
tion changes depending on the choice of the discretization, the phase transition 
might be seen only in the lattice models, and as a consequence we can obtain 
no information on the phase transition of the continuous model by using the 
numerical simulations. Hence, it is natural to consider that the order of the 
transition depends on the discretization of Hamiltonian. Therefore, it is worth- 
while to investigate further the tethered model defined by Hamiltonian which 
is different from the above-mentioned bending energy and Gaussian term. 

In this paper, we studied three types of tethered models on triangulated spher- 
ical surfaces, which are obtained by dividing the icosahedron as in [26]. The 
Hamiltonian of the first model contains the ordinary Gaussian energy and a 
bending energy, which is defined by the normal vector of vertices. Hence, the 
discretization of the bending energy in the first model is slightly different from 
that of the standard bending energy. It was previously reported that the first 
model undergoes a discontinuous transition on Voronoi lattices [25] . However, 
it remains unclear whether the discontinuous transition can be seen on any 
other lattices. Thus, we can also check that the discontinuous transition is 
independent of the lattice structure by investigating the first model defined 
on lattices constructed from the icosahedron. 

The second model is called the Nambu-Goto surface model [19,21,31]. At 



2 



present, it is unclear whether the Nambu-Goto model undergoes a discon- 
tinuous phase transition when the model is defined on surfaces with extrinsic 
curvature. The area energy, which is given by the sum of the area of the tri- 
angles, corresponds to the Gaussian energy in the first (HPK) model. It is 
well-known that the Nambu-Goto model without the bending energy term is 
not well-defined due to the lack of the smooth phase and the appearance of 
spikes on the surface [36]. The area term imposes a constraint only on the 
area of the triangles, and hence allows the appearance of spiky configurations 
made up of oblong triangles. On the contrary, some additional terms in the 
Hamiltonian can make the model well-defined. In fact, it was reported recently 
that the Nambu-Goto model with a deficit angle term, which is an intrinsic 
curvature, is well-defined and undergoes a discontinuous transition between 
the smooth phase and a tubular phase [21]. The second model in this pa- 
per is also well-defined [31,19], because the Hamiltonian includes a bending 
energy, which is an extrinsic curvature defined according to the dual lattice 
formulation of the discrete mechanics by Lee [37] . 

The third model is a tensionless model whose Hamiltonian contains the stan- 
dard bending energy and a hard-wall potential which gives an upper bound on 
the bond length. There is no Gaussian term in the Hamiltonian. The bending 
energy of the model is identical to the one in [26] . Therefore, the discretization 
of the third model is the same as in [26]. However, it is not yet clear whether 
a tensionless tethered model undergoes a discontinuous transition. 

It is nontrivial whether those three discrete models lead to the same result on 
the order of the transition. The standard bending energy for the HPK model 
is not equivalent to the above mentioned bending energy for the Nambu-Goto 
surface model. In fact, the Nambu-Goto model with the standard bending en- 
ergy is ill-defined because of the lack of the smooth phase and the appearance 
of spikes in the whole region of the bending rigidity [36]. Thus, it is worthwhile 
to study numerically whether these three models have the same continuous 
limit for the order of the transition. 



We will see that the three models undergo a discontinuous transition on spher- 
ical surfaces relatively larger than those used to date, and then we will un- 
derstand that the discontinuous transition of the model is independent of the 
discretization of the Hamiltonian. Therefore, the results obtained in this pa- 
per together with those in [26] lead us to conclude that a first-order phase 
transition can be observed in a spherical tethered surface model defined by 
Hamiltonian that includes a bending energy term. This means that the first- 
order transition of triangulated surfaces is universal in the sense that the order 
of the transition is independent of discretization of the Hamiltonian. 
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2 Models 



The partition function of the models is defined by 

N 



[ fldXiexp[-SiX)], (1) 

i=i 



where is the total number of vertices, S{X) is Hamiltonian of the model. 
The center of the surface is fixed to remove the translational zero-mode. The 
self- avoiding interaction is not assumed in the Hamiltonian, hence the surface 
is considered to be phantom. 

We investigate three types of model which differ from each other in energy 
terms included in the Hamiltonian S{X). The first model denoted by model 1 
is defined by 



S{X)^Si + hS2, 



N 



S, = J2{X,-Xjf, 52 = EE[l-n(^)-^.»]' (model 1) (2) 

(id) *=i 

where b is the bending rigidity, nj(j) is the unit normal of the triangle 
meeting at the vertex i. The symbol n(i) in Eq. (2) is a normal vector of the 
vertex i, and it is defined by 



N- 



where J2j{i) denotes the summation over triangles linked to the vertex i, 
and ^zi^(j) is the area of Model 1 is considered a HPK model as mentioned 
in the introduction. 

The second model denoted by model 2 is defined by 
S{X) ^Si + bS2, 

S^ = J2A^, 52 = E^E^'(l-cos^.), (model 2) (4) 

A A i=l 

where A/^ is the area of the triangle A, 1^ is the edge length of A, and 9i is 
the angle between the triangles attaching to the edge Zj. 

Figures 1(a), (b) show the triangles whose normal vectors are interacting with 
one of the shaded triangles corresponding to model 1 and model 2, respectively. 
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(a) (b) 

Fig. 1. The range of interaction of the shaded triangle in (a) model 1, and (b) model 
2. in Eq. (4) is the area of the shaded traiangle in (b), and li is the edge length 
of the triangle. 

in Eq. (4) is the area of the shaded triangle in Fig. 1(b), and li in Eq. (4) 
is the edge length of the triangle. 

It should be noted that the Nambu-Goto model is considered a classical ana- 
logue of the Polyakov string model. In fact, the Polyakov string model reduces 
to the Nambu-Goto model if the metric, which is a dynamical variable of the 
Polyakov string, is fixed to the induced metric [38]. 

The third model denoted by model 3 is defined by the Hamiltonian 

Six) = bS2 + V, 

52 = 5:(l-n,-n,), V ^Y.^{\Xi-X,\), (models). (5) 

S2 in S{X) of Eq. (5) is the bending energy, which was already introduced in 
the Introduction, and is a hard-wall potential, which gives an upper bound 
on the bond length. Zl(ij) in V denotes the sum over bonds (i, j) connecting 
the vertices i and j, and Y.^j in 'S'2 is the sum over triangles i,j sharing a 
common bond. The symbol n, in S'2 denotes the unit normal vector of the 
triangle i. The range of interaction of Hj in S'2 is identical to that of model 
2, shown in Fig. 1(b). The symbol V(|Xj — Xj|) in Eq. (5) is the potential 
between the vertices i and j, and is defined by 

V{\X,-X,\) = 



The value of tq in the right hand side of Eq. (6) is fixed to tq = vTT. Then 
we have (^{Xi — XjY) / N — 3/2, which is satisfied when the Gaussian term 
5*1 = Y.iX—XjY is included in the Hamiltonian without the hard- wall potential 




(0<|X,-X,|<ro), 

(6) 

00 (otherwise) . 
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V. 



If it were not for the constraint — Xj| < ro, the size of the surface grows 
larger and larger in the MC simulations. Thus the constraint \Xi—Xj \ <ro is 
necessary for the bond length to have a well-defined value when the Gaussian 
term 5*1 is not included in the Hamiltonian. We note that the model seems to 
depend on a hidden length scale introduced by Tq, and as a consequence the 
model appears to be ill-defined. However, we checked in Ref. [20] that there is 
no ro-dependence on the results. This is a consequence of the scale invariant 
property of the model, which will be mentioned in the next section. Therefore, 
we use ro = vTT in the MC simulations. 

We should comment on a relationship between the discrete bending energy 
5*2 in Eqs. (2), (4), and (5) and the continuous one of the Polyakov-Kleinert 
string model. The functional action of the Polyakov-Kleinert string in is 
defined by 

S(X,g) = I / + J ^sdhKlKt, (7) 

M M 

where {gab) is the first fundamental form (or the metric tensor) of the surface 
M. The K^^ is defined by K\ = g^'^Kac, where Kab is the second fundamental 
form (or the extrinsic curvature tensor) on M in [5]. In S{X, g), X denotes 
a mapping from the surface M to R^, where the image X[M) is the membrane. 

It is easy to sec that (1/2) cfxK^K^ in the second term of S{X, g) in Eq. 

(7) is equivalent to 5*2 in Eqs. (2), (4), and (5) if the metric tensor [gab) is given 
by the induced one; gab = ^f^^r- In fact, the second fundamental form K^b 
can be expressed by Kab = n^dadbX'^ = —dafi^dbX^, where is the unit nor- 
mal vector of the surface. Then by using this relation and the induced metric, 
the second term in S{X,g) can be expressed by (1/2) ^Jg dP'xg"''^dan^di,n^ . 
This action is identical to the one of the cr-model, and hence a natural dis- 
cretization of this term is given by S2 in Eqs. (2), (4), and (5). 



3 Monte Carlo technique 

The canonical Monte Carlo (MC) technique is used to update the variable X. 
The new position X'- of the vertex i is given by X'- — Xi-\-AX, where AX is 
chosen randomly in a small sphere for model 1 and model 2. The radius of the 
small sphere is chosen at the start of the MC simulations to maintain about 
50% acceptance rate. The new position X[ is accepted with the probability 
Min[l,exp(-Z\S')], where AS is given by Zl5' = 5'(new) -^(old). 
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Under the constraint of Eq.(6) for the edge length for model 3, the new position 
X- is accepted with the probability Min[l, exp {—AS)]; namely the Metropolis 
accept /reject procedure is applied only for satisfying the constraint of 
Eq.(6). This constraint for the edge length is imposed by Tq = 1.1 as described 
in the previous section. The rate of acceptance for the constraint of Eq. (6) is 
about 50%. 

A random number called Mersenne Twister [39] is used in the MC simulations. 
We use two-sequences of random numbers; one for the 3-dimensional move of 
vertices X and the other for the Metropolis accept/reject in the update of X. 

For model 1 and model 2, the minimum bond length is not assumed, while 
the minimum area of triangle is assumed to be 10~^ x ^40, where is the 
mean area of triangles computed at every 250 MCS (Monte Carlo sweeps) 
and is almost constant throughout the MC simulations. The area of almost all 
triangles generated in the simulations is larger than the lower bound 10""^ xAq. 

For model 3, the minimum area of triangle is assumed to be 1 x 10~^. The 
area of almost all triangles generated in the MC simulations is larger than the 
lower bound 1 x 10~^. The minimum bond length is assumed to be 1 x 10~^, 
however, we find no bond whose length is less than the lower bound 1 x 10~^ 
in the MC simulations. This means that no constraint is imposed on the lower 
bound of the bond length. 

The Hamiltonians in Eqs. (2), (4), and (5) are defined on the surfaces which 
are constructed uniformly by the co-ordination number. They are obtained 
by dividing the icosahedron. By dividing every edge of the icosahedron into 
L-pieces of the same length, we have a triangulated lattice of size iV = lOL^+2. 
These lattices are characterized by iV5 = 12 and A''6 = A^ — 12, where Ng is the 
total number of vertices with co-ordination number q; thus we have lattices in 
which 12 vertices are of qi — 5, and all other vertices Qi — Q. 

We comment on the unit of physical quantities in the model. The scale of 
length in the model can be arbitrarily chosen because of the scale invariant 
property of the partition function in Eq. (1). Then, by letting a be a length 
unit in the model, we can express all quantities with unit of length in terms 
of a. Hence, the unit of 5*1 is a^. Let A be the surface tension coefficient, then 
S in Eqs. (2), (4) for model 1, model 2 can be written as S = XSi + bS2. 

Thus, the unit of A can be written as kT/a^, where k is the Boltzmann con- 
stant, T is the temperature. The unit of h is therefore expressed by kT. The 

coefficient A is assumed to be A = 1 in this paper. This is always possible 
because of the scale invariant property of the partition function of the model. 
The unit of A can be influenced by the length unit a of the model, and the 
value of A varies depending on a. 
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However, A can always be fixed to A= 1. This is the reason why the length unit 
a is arbitrary chosen in the model. Note also that varying the temperature T 
is efi:ectively identical to varying h in the model. The dependence of A on T 
can be absorbed into a redefinition of A, which can be fixed to A = 1 due to the 
scale invariant property as stated above. The parameter ro in Eq. (6) plays a 
role in the length unit a of model 3. 



4 Results 

First, we summarize the total number of MCS iterated in the simulations. At 
the transition point of model 1, 1.8x10^, 1.8x10^, and 1.5x10^ MCS were 
iterated on the surfaces of = 6252, N = 4412, and N = 2562, respectively. 
At the transition point of model 2, 7.8 x 10^ 7x 10^, and 6x 10^ MCS were 
iterated on the surfaces of A^= 10242, A^ = 6762, and A" = 4842, respectively. 
For model 3, the total number of MCS at the transition point was about 3x10^, 
5xl0^ and 6.5x10^ for the surfaces of A^ = 4842, A^ = 8412, and A^ = 16812, 
respectively. A relatively small number of MCS were performed at b which is 
not the transition point on each surface. 

The mean square size is defined by 

i i 

where X is the center of the surface. 

Figure 2(a) is a normalized distribution h{X^) (or the histogram) of X"^ for 
model 1 on the surface of size A^ = 4412 at the transition point 6 = 0.4565. The 
variation of against MCS is plotted in Fig. 2(b). Figures 2(c) and 2(d) are 
those obtained on the surface of A^ = 6252 at 6 = 0.455. Two distinct peaks 
visible in h{X'^) imply that the size of surfaces in one phase is very different 
from that in the other phase at the transition point. 

Figure 3(a) is a normalized histogram h{X'^) of X^ for model 2 on the surface 
of size A^ = 6762 at the transition point 6 = 0.1144. The variation of against 
MCS is plotted in Fig. 3(b). Figures 3(c) and 3(d) are those obtained on the 
surface of A^ = 10242 at 6=0.1144. 

Figure 4(a) is a normalized histogram h{X'^) of for model 3 on the surface 
of size X = 8412 at the transition point 6 = 0.686. The variation of against 
MCS is plotted in Fig. 4(b). Figures 4(c) and 4(d) are those obtained on the 
surface of A^= 16812 at 6=0.683. 
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Fig. 2. (a) The normalized histogram h{X^) of for model 1 on the N = 4412 
surface at 6=0.4565, and (b) the variation of X'^ against MCS. (c) h{X^), and (d) 
the variation of X^ obtained on the N = 6252 surface at 6 = 0.455. Two distinct 
peaks on each h{X'^) imply that the size of surfaces in one phase is very different 
from that in the other phase at the transition point of model 1. 




400^^2 6Q0 
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Fig. 3. (a) The normalized histogram h{X'^) of X^ for model 2 on the N = 6762 
surface at 6 = 0.1144, and (b) the variation of against MCS, (c) h{X^) and (d) 
the variation of X"^ obtained on the A'" = 10242 surface at 6 = 0.1144. 

In order to obtain the Hausdorff dimensions i^smo (-f^cru) in the smooth (crum- 
pled) phase close to the transition point, the mean value of was obtained 
by averaging X"^ over a small region at each peak of h{X'^). For model 1, we 
assumed the region: 4 < < 10 and 14 < < 19 at 6= 0.458 on the N = 1442 
surface, 6 < < 18 and 20 < < 31 at 6 = 0.456 on the = 2562 surface, 
7 < X2 < 25 and 35 < < 52 at 6 = 0.4565 on the N = 4142 surface, and 
8<X2<30 and 50<X2<72 at 6 = 0.455 on the TV = 6252 surface. 



For model 2, we assumed the region: 40 < X^ < 90 and 100 < < 150 at 
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(c) N=16812 
b=0.683 
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Fig. 4. (a) The normalized histogram h{X^) of X'^ for model 3 on the N = 8412 
surface at 6 = 0.686, and (b) the variation of X'^ against MCS, (c) h{X'^) and (d) 
the variation of X"^ obtained on the iV = 16812 surface at 6 = 0.683. 

b = 0.115 on the N = 2562 surface, 50 < < 150 and 170 < < 270 at 
b = 0.1145 on the = 4842 surface, 65 < < 220 and 230 < < 380 at 
6 = 0.1144 on the iV = 6762 surface, and 120<X2<320 and 360<X2<560 at 
6 = 0.1144 on the iV= 10242 surface. 



For model 3, we assumed the region: 8<X^<21 and 23<X^<36 at 6 = 0.692 
on the surface of iV = 2562, 9 < < 35 and 39 < < 64 at 6 = 0.688 on 
TV = 4842, 11 < X2 < 49 and 70 < X^ < 108 at 6 = 0.686 on X = 8412, and 
16<X2<55 and 135<X2< 205 at 6 = 0.683 on the surface of X= 16812. 




(a) N (b) N (c) N 



Fig. 5. Log-log plots of X^ against N obtained in the smooth phase and in the 
crumpled phase close to the transition point for (a) model 1, (b) model 2, and (c) 
model 3. The error bars on the data represent the standard deviations. 

Figures 5(a),(b),(c) are log-log plots of X^ against X for model 1, model 2, 
and model 3 respectively, where the values of X^ are obtained by averaging X^ 
in the regions described above. Error bars on the data represent the standard 
deviations. The straight lines are drawn by fitting the data X^ to 

X^ ~ X^/^. (9) 
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Then we have 




2.22 ±0.13 (modell), 

2.11 ±0.18, //cru(2) = 2.39 ± 0.51 



(model 2), 



(10) 



^smo(3) = 2.13 ± 0.11 (model 3). 

where i in i?smo,cru(^) denotes the model i{— 1,2,3), and smo (cm) denotes 
the smooth (crumpled) phase, ifsmo(l), ^^smo(2), and i?smo(3) in Eq. (10) are 
comparable to the topological dimension H = 2 within the errors. This data 
indicates that the surfaces are smooth in the smooth phase at the transition 
point for the three models. 

-f^cru(2) in Eq.(lO) is less than the physical bound H — 3. This implies that 

the surface is not completely crumpled but relatively smooth in the crumpled 
phase at the transition point of model 2. A similar phenomenon occurs for the 
standard tethered model [26]. Note also that i/cru(2) in Eq. (10) is comparable 
to the theoretical prediction H — 2.39(23) within the error, which corresponds 
to the scaling exponent i/ = 0.84±0.04 [40] where v — 2/H. 

On the contrary, we got ifcru(l) = 3.27±0.88 in the crumpled phase at the 
transition point of model 1 by fitting the data to Eq. (9). However, it seems 
that ifcru(l) is not well-defined because of the large error. The logarithmic 
divergence [3,41,42,43]; = co + CilogA^, is also expected for X'^ in the 
crumpled phase. In fact, the residual sum of squares RSS is RSS = 0.424 
for the log-log fit of X"^ in the crumpled phase at the transition point, whereas 
the linear-log fit gives RSS = 0.0906 for these X^, where RSS is defined by 
RSS = X^[(data — fitting formula) /error] ^. The fact that the hnear-log fit is 
better than the log-log fit indicates that scales according to the logarithmic 
divergence and that i7cru(l) = 3.27±0.88 is not well-defined. 

We also have i7cru(3) =4.53± 1.23 by fitting X^ in the crumpled phase at the 
transition point of model 3. The dashed line in Fig. 5(c) is drawn by using 
this i7ciu(3), and indicates that ifcru(3) is ill-defined for the same reason as for 
H„^{1). In fact, we have RSS = 0.00438 for the linear-log fit and RSS = 0.221 
for the log- log fit in Fig. 5(c). Thus, the linear- log fit is better than the log- log 
fit for X^ in the crumpled phase of model 3. This allowed us to conclude that 
-f^cru(3) = 4.53 ± 1.23 is an ill-defined value. 

We understand also that a choice of discretization of the Hamiltonian infiu- 
ences the scaling property of X^ in the crumpled phase at the transition point. 

Figure 6(a) is a normalized histogram h{S2) of S2/NBB obtained on the surface 
of size X = 4412 at the transition point 6 = 0.4565, where Nbb denotes Nbb = 
2Nb, and Nb is the total number of bonds given by Nb = 3X — 6. The 
bending energy 5*2 of model 1 is divided by Nbb so that S2/NBB can be 
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Fig. 6. (a) The normalized histogram h{S2) of S2/NBB, and (b) the variation of 
S2/NBB against MCS, on the iV = 4412 surface at 6 = 0.4565, (c) h{S2) and (d) the 
variation of S2/NBB on the AT = 6252 surface at 6=0.455. Nbb = 2Nb. 

regarded as 1 — (cos^). Two peaks can be seen clearly in h{S2). The variation 

of 5*2 against MCS is plotted in Fig. 6(b). Figures 6(c) and 6(d) show h{S2) 
and the variations of S2 obtained on the surface of iV = 6252 at 6 = 0.455. 
We also find two distinct peaks in h{S2) in Figs. 6(c). These two peaks in 
Figs. 6(a) and 6(c) obviously show that model 1 undergoes a first-order phase 
transition. Note also that the discontinuous transition is independent of the 
lattice structure. It was reported in [25] that the same model undergoes a 
discontinuous transition on Voronoi triangulated spherical surfaces, which are 
different from those used in this paper. 




Fig. 7. (a) h{S2) of S2/NT, and (b) the variation of S2/NT against MCS, on the 
N = 6762 surface at 6 = 0.1144, (c) h{S2) and (d) the variation of S2/NT on the 
A^ = 10242 surface at 6=0.1144. Nt = 2N-4. 

Figure 7(a) is a normalized histogram h{S2) of S2/NT obtained on the A'" = 6762 
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surface at the transition point 6 = 0.1144, where A^^ denotes the total number 
of triangles given by Nt = 2N—4. The variation of 5*2 against MCS is plotted in 
Fig. 7(b). Figures 7(c) and 7(d) show h{S2) and the variations of 5*2 obtained 
on the = 10242 surface at 6 = 0.1144. The double peak structures of h{S2) in 
Figs. 7(a) and 7(c) obviously show that model 2 undergoes a first-order phase 
transition. 




Fig. 8. (a) h{S2) of S2/NB, and (b) the variation of S2/NB against MCS, on the 
N = 8412 surface at b = 0.686, (c) h{S2) and (d) the variation of S2/NT on the 
iV = 16812 surface at 6 = 0.683. Nb = 3N-6. 

Figure 8(a) is a normalized histogram h{S2) of S2/NB obtained on the N = 
8412 surface at the transition point 6 = 0.686. The variation of 5*2 against MCS 
is plotted in Fig. 8(b). Figures 8(c) and 8(d) show h{S2) and the variations of 
5*2 obtained on the A^ = 16822 surface at 6 = 0.683. The double peak structures 
in h{S2) in Figs. 8(a) and 8(c) obviously show that model 3 undergoes a first- 
order phase transition. 

Figures 9(a), (b), (c) show X"^ against 6 obtained in model 1, model 2, and 
model 3, respectively. Figure 9(d) shows the bending energy S2/NBB against 
6 obtained in model 1. Figure 9(e) shows X"^ and S2/NT against 6 obtained 
in model 2, and Fig. 9(f) shows X"^ and S2/NB against 6 obtained in model 
3. Note that the value for 5*2 of model 2 is very different from that for model 
1 or model 3. The reason for this is because S'2/A^T can not be regarded as 
1 — (cos^), which can be demonstrated in Eq. (4). 

X^, S2/NBB, and S2/NT in Figs. 9(a), (b), (d), (e) are almost smooth, and no 
discontinuous change is seen for intermediate bending rigidity 6. The reason 
for those smooth behaviors is because the surface size is relatively small for 
model 1 and model 2. We expect that the discontinuous changes can be seen 
in those quantities obtained on large surfaces. However, we consider that the 
sizes N = 6252 and N = 10242 are sufficiently large, because the discontinuous 
nature of the transition is obvious from the double peak in h{S2) in Figs. 6 
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Fig. 9. The mean square size against the bending rigidity b of (a) model 1, (b) 
model 2, (c) model 3, and the bending energy (d) S2/NBB against h of model 1, (e) 
S2/NT against b of model 2, (f) S2/NB against b of model 3. Nbb is twice that of 
Nb the total number of bonds and Nt is the total number of triangles. 

and 7. On the contrary, we see discontinuous changes in X^, S2/NB shown in 
Figs. 9(c), (f) for model 3. Surface size A^= 16812 seems sufficiently large for 
model 3 showing the discontinuity of X^, S2/NB. 
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Fig. 10. The specific heat C52 against the bending rigidity 6 of (a) model 1, (b) 
model 2, and (c) model 3. Log-log plot of the peak values C^^^ against N of (d) 
model 1, (e) model 2, and (f) model 1. The straight lines in (d) and (e) are drawn 
by fitting the largest three Cg^^^ to Eq. (12). The solid (dashed) line in (f) is drawn 
by fitting the largest (smallest) four to Eq. (12). 
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The specific heat is defined by 

C& = ^((^2-(^2))'). 



(11) 



Figure 10(a) shows of model 1 on the surfaces of size A'" = 6252, N = AAV2, 
and A^ = 2562. Sharp peaks in Fig. 10(a) indicate that model 1 undergoes a 
discontinuous transition. Figure 10(b) shows of model 2 on the surfaces 
of iV= 10242, A^ = 6762, and iV = 2562, and Fig. 10(c) shows Cs^ of model 3 
on the surfaces of 16812, Ar = 8412, A^ = 4842, and A^ = 2562. 

Figure 10(d) is a log-log plot of the peak value C^g^^ against including the 
results obtained on the = 1442 surface for model 1. The peaks C^^^^ are 
plotted in Fig. 10(e) and in Fig. 10(f) for model 2 and model 3, respectively. 

The straight lines in Figs. 10(d), (e) were drawn by fitting the largest three 

c^r to 

~ N% (12) 



where ex is a critical exponent of the transition. The straight line denoted by 
large four {small four) in Fig. 10(f) were drawn by fitting the largest (smallest) 
four C^^"" to Eq. (12). 

Thus, we have 



(7i = 0.93 ± 0.28 (model 1), 

(72 = 0.83 ±0.13 (model 2), (13) 
(73 = 1.008 ± 0.048, (model 3), 

and (73(small) = 0.45 ±0.14 [N < 2562) for model 3. The exponents (7i, (72, (73 
in Eq. (13) indicates that the phase transition is of the first order in all three 
models. On the contrary, (73(small) = 0.45(14) for [N < 2562) indicates that 
model 3 appears to undergo a continuous transition on small-sized surfaces 
of iV < 2562. The result a^{sm.al\) = 0.45(14) is almost comparable to cr = 
0.406(14) of the model on Voronoi triangulated spherical surfaces of size up 
to N — 1500 in [20]. Two different behaviors for against N shown in Fig. 
10(f) are consistent with the fact that two distinct peaks in h{S2) are observed 
only on large (A/">4842) surfaces. 

Figures 11(a), (b),(c) are the autocorrelation coefficient A{X'^) of X'^ defined 

by 

A{X') = ^'f^^^f^^'\ -n + nx 500, (n = 1, 2, . . .), (14) 
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Fig. 11. Autocorrelation coefficient A{X'^) of obtained on (a) the model 1 surface 
of Ar = 6762 at 6 = 0.454(A), 6 = 0.455(0), and 6 = 0.457(n), (b) the model 2 surface 
of iV = 10242 at 6 = 0. 1138(A), 6 = 0.1144(0), and 6 = 0. 1148(0), and (c) the model 
3 surface iV = 16812 at 6 = 0.682(A), 6=0.683(0), and 6 = 0.684(0). 

from which we can see the convergence speed of in the simulations. A{X'^) 
in Fig. 11(a) are obtained on the model 1 surface of = 6762 at 6 = 0.454, 
6 = 0.455, and 6 = 0.457, and those in Fig. 11(b) are obtained on the model 
2 surface of = 10242 at 6 = 0.1138, 6 = 0.1144, and 6 = 0.1148. Figure 
11(c) shows A(X2) obtained on the model 3 surface A^ = 16812 at 6 = 0.682, 
6 = 0.683, and 6 = 0.684. In these figures, wc find the phenomenon of critical 
slowing down, which reflects the transition. 

The horizontal axes in Figs. 11(a), (b),(c) represent 500 x n {n = 1,2,---)- 
MCS, which is a sampling-sweep between the samples X^(rj) and X^(rj+i). 
The critical slowing down is clearly seen in the figures. The reason for this 
critical slowing down is because the volume of phase space (C R-^) for X at 
the transition point becomes larger than those at the crumpled phase and at 
the smooth phase. It is remarkable that the phase space volume in the smooth 
phase seems comparable to that in the crumpled phase, despite the fact that 
the surfaces in the smooth phase extend to large regions in R3. This can also 
be seen in a tethered model with intrinsic curvature on a sphere [33] and in a 
tethered model on a torus [23] . 

Figure 12(a) is a snapshot of model 1 surface of = 6762 in the crumpled 
phase at 6 = 0.455, and Fig. 12(b) is the one in the smooth phase at the same 
6. The sections of the surface shown in Figs. 12(a) and 12(b) are depicted in 
Figs. 12(c) and 12(d), respectively. Snapshots were drawn at the same scale. 

Figure 13(a) is a snapshot of model 2 surface of A^= 10242 in the crumpled 
phase at 6=0.1144, and Fig. 13(b) is the one in the smooth phase at the same 
6. The sections of the surface shown in Figs. 13(a) and 13(b) are depicted in 
Figs. 13(c) and 13(d), respectively. Snapshots were drawn at the same scale. 

Figure 14(a) is a snapshot of the iV = 16812 surface in the crumpled phase at 
6 = 0.683, and Fig. 8(b) is the one in the smooth phase at the same 6. The 
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Fig. 12. Snapshots of model 1 surfaces at (a) the crumpled phase and at (b) the 
smooth phase, and (c) the section of the surface in (a), and (d) the section of the 
surface in (b). The snapshots were obtained at the transition point 6 = 0.455 on the 
surface of A^ = 6762. 



mean square size is about X"^ = 39 and = 178 in Figs. 14(a) and 14(b), 
respectively. The sections of the surface in Figs. 14(a) and 14(b) are depicted 
in Figs. 14(c) and 14(d), respectively. 

Surfaces are rough at short scales even for the smooth phase shown in Figs. 
12(b), 13(b), 14(b), whereas they are smooth at the long range scale. Surfaces 
rough at short scales can also be seen deep in the smooth phase. Only a 
spherical monolayer surface is apparent in the smooth phase. None of non- 
spherical surfaces; such as oblong surfaces, linear ones, and branched-polymer 
like ones, were observed. The surfaces in the crumpled phase at the transition 
point are completely collapsed in both model 1 and model 3. These crumpled 
surfaces are in contrast to those in model 2 and those in [26], where the 
crumpled state at the transition point seems not completely crumpled and 
seems to be characterized by a Hausdorff dimension less than the physical 
bound H = 3. 
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Fig. 13. Snapshots of model 2 surfaces at (a) the crumpled phase and at (b) the 
smooth phase, and (c) the section of the surface in (a), and (d) the section of the 
surface in (b). The snapshots were obtained at the transition point 6 = 0.1144 on 
the surface of A^ = 10242. 

5 Summary and conclusions 



This paper aimed to show that the first-order transition of the tethered surface 
model for biological membranes is independent of the discretization of the 
Hamiltonian. Three types of tethered surface models were investigated. The 
first model (model 1) is a discrete model of Helfrich and Polyakov-Kleinert, 
the second (model 2) is that of Nambu-Goto, and the third (model 3) is a 
tensionless model of Helfrich and Polyakov-Kleinert. The first and the second 
models are closely related to each other in the context of the string model. 



The canonical MC simulations were carried out on spherical surfaces of size 
up to iV = 6252 for model 1, A^ = 10242 for model 2, and A^= 16812 for model 
3, where the lattices were constructed by dividing the icosahedron. A discrete 
form of the Hamiltonian of model 1 is a linear combination of the Gaussian 
term and a bending energy term, which is defined by using a normal vector 
of vertices. The Hamiltonian of model 2 contains an area term and a bending 
energy term, which makes the model well-defined. The Hamiltonian of model 
3 contains a bending energy term and a hard-wall potential, which gives an 
upper bound on the bond length. The first-order transition of models 1, 2, 
and 3 were observed on lattices of size N > 2562, > 6762, and > 4842, 
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Fig. 14. Snapshots of model 3 surfaces at (a) the crumpled phase and at (b) the 
smooth phase, and (c) the section of the surface in (a), and (d) the section of the 
surface in (b). The snapshots were obtained at the transition point 6 = 0.683 on the 
surface of iV = 16812. 

respectively. 

We have shown in this paper that the first-order crumpling transition is uni- 
versal on spherical surfaces of the extrinsic curvature models for biological 
membranes. In fact, we have found that the first-order transition is indepen- 
dent of a choice of Hamiltonian of the tethered surface model, if it is discretized 
on a sphere. Therefore, the results obtained in this paper together with those 
in [26] lead us to conclude that a first-order phase transition can be observed 
in a spherical tethered surface model whose Hamiltonian includes a bending 
energy term. It is also found that the Hausdorff dimension H in the smooth 
phase at the transition point is almost identical to the topological dimension 
H = 2 and is also independent of the discretization of the Hamiltonian. On the 
contrary, a choice of discretization of the Hamiltonian influences the scaling 
property of in the crumpled phase at the transition point. Nevertheless, 
we consider that the crumpled state of the model is not a phantom one but 
a real physical one, because the Hausdorff dimension of the crumpled phase 
in model 2 is iJcru(2) = 2.39(51), which is less than the physical bound and is 
also comparable to the result in [26]. 

Finally, we emphasize that our results suggest that biological membranes can 
exhibit the crumpling transition. Experimental investigations on the crum- 
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pling transition are expected. Further numerical studies on the transition 
should also be performed on larger surfaces. 

This work is supported in part by a Grant-in-Aid of Scientific Research, No. 
15560160. 
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